Rare variant analyses in large-scale cohorts identified SLC13A1 associated with chronic pain.

ABSTRACT: Chronic pain is a prevalent disease with increasing clinical challenges. Genome-wide association studies in chronic pain patients have identified hundreds of common pathogenic variants, yet they only explained a portion of individual variance of chronic pain. With the advances in next-generation sequencing technologies, it is now feasible to conduct rarer variants studies in large-scale databases. Here, we performed gene-based rare variant analyses in 200,000 human subjects in the UK biobank whole-exome sequencing database for investigating 9 different chronic pain states and validated our findings in 3 other large-scale databases. Our analyses identified the SLC13A1 gene coding for sodium/sulfate symporter associated with chronic back pain and multisite pain at the genome-wide level and with chronic headache, knee, and neck and shoulder pain at the nominal level. Seven loss-of-function rare variants were identified within the gene locus potentially contributing to the development of chronic pain, with 2 of them individually associated with back pain and multisite pain. These 2 rare variants were then tested for replication in 3 other biobanks, and the strongest evidence was found for rs28364172 as an individual contributor. Transcriptional analyses of Slc13a1 in rodents showed substantial regulation of its expression in the dorsal root ganglia and the sciatic nerve in neuropathic pain assays. Our results stress the importance of the SLC13A1 gene in sulfate homeostasis in the nervous system and its critical role in preventing pain states, thus suggesting new therapeutic approaches for treating chronic pain in a personalized manner, especially in people with mutations in the SLC13A1 gene.

Implementation of a graph-embedded topic model for analysis of population-level electronic health records.

To address the need for systematic investigation of the phenome enabled by ever-growing genotype and phenotype data, we describe our step-by-step software implementation of a graph-embedded topic model, including data preprocessing, graph learning, topic inference, and phenotype prediction. As a demonstration, we use simulated data that mimic the UK Biobank data as in our original study. We will demonstrate topic analysis to discover disease comorbidities and computational phenotyping via the inferred topic mixture for each subject. For complete details on the use and execution of this protocol, please refer to Wang et al. (2022).1.

Acute inflammatory response via neutrophil activation protects against the development of chronic pain.

The transition from acute to chronic pain is critically important but not well understood. Here, we investigated the pathophysiological mechanisms underlying the transition from acute to chronic low back pain (LBP) and performed transcriptome-wide analysis in peripheral immune cells of 98 participants with acute LBP, followed for 3 months. Transcriptomic changes were compared between patients whose LBP was resolved at 3 months with those whose LBP persisted. We found thousands of dynamic transcriptional changes over 3 months in LBP participants with resolved pain but none in those with persistent pain. Transient neutrophil-driven up-regulation of inflammatory responses was protective against the transition to chronic pain. In mouse pain assays, early treatment with a steroid or nonsteroidal anti-inflammatory drug (NSAID) also led to prolonged pain despite being analgesic in the short term; such a prolongation was not observed with other analgesics. Depletion of neutrophils delayed resolution of pain in mice, whereas peripheral injection of neutrophils themselves, or S100A8/A9 proteins normally released by neutrophils, prevented the development of long-lasting pain induced by an anti-inflammatory drug. Analysis of pain trajectories of human subjects reporting acute back pain in the UK Biobank identified elevated risk of pain persistence for subjects taking NSAIDs. Thus, despite analgesic efficacy at early time points, the management of acute inflammation may be counterproductive for long-term outcomes of LBP sufferers.

A graph-embedded topic model enables characterization of diverse pain phenotypes among UK biobank individuals.

Large biobank repositories of clinical conditions and medications data open opportunities to investigate the phenotypic disease network. We present a graph embedded topic model (GETM). We integrate existing biomedical knowledge graph information in the form of pre-trained graph embedding into the embedded topic model. Via a variational autoencoder framework, we infer patient phenotypic mixture by modeling multi-modal discrete patient medical records. We applied GETM to UK Biobank (UKB) self-reported clinical phenotype data, which contains 443 self-reported medical conditions and 802 medications for 457,461 individuals. Compared to existing methods, GETM demonstrates good imputation performance. With a more focused application on characterizing pain phenotypes, we observe that GETM-inferred phenotypes not only accurately predict the status of chronic musculoskeletal (CMK) pain but also reveal known pain-related topics. Intriguingly, medications and conditions in the cardiovascular category are enriched among the most predictive topics of chronic pain.

Detecting multi-way epistasis in family-based association studies.

The era of genome-wide association studies (GWAS) has led to the discovery of numerous genetic variants associated with disease. Better understanding of whether these or other variants interact leading to differential risk compared with individual marker effects will increase our understanding of the genetic architecture of disease, which may be investigated using the family-based study design. We present M-TDT (the multi-locus transmission disequilibrium test), a tool for detecting family-based multi-locus multi-allelic effects for qualitative or quantitative traits, extended from the original transmission disequilibrium test (TDT). Tests to handle the comparison between additive and epistatic models, lack of independence between markers and multiple offspring are described. Performance of M-TDT is compared with a multifactor dimensionality reduction (MDR) approach designed for investigating families in the hypothesis-free genome-wide setting (the multifactor dimensionality reduction pedigree disequilibrium test, MDR-PDT). Other methods derived from the TDT or MDR to investigate genetic interaction in the family-based design are also discussed. The case of three independent biallelic loci is illustrated using simulations for one- to three-locus alternative hypotheses. M-TDT identified joint-locus effects and distinguished effectively between additive and epistatic models. We showed a practical example of M-TDT based on three genes already known to be implicated in malaria susceptibility. Our findings demonstrate the value of M-TDT in a hypothesis-driven context to test for multi-way epistasis underlying common disease etiology, whereas MDR-PDT-based methods are more appropriate in a hypothesis-free genome-wide setting.

Human TYK2 deficiency: Mycobacterial and viral infections without hyper-IgE syndrome.

Autosomal recessive, complete TYK2 deficiency was previously described in a patient (P1) with intracellular bacterial and viral infections and features of hyper-IgE syndrome (HIES), including atopic dermatitis, high serum IgE levels, and staphylococcal abscesses. We identified seven other TYK2-deficient patients from five families and four different ethnic groups. These patients were homozygous for one of five null mutations, different from that seen in P1. They displayed mycobacterial and/or viral infections, but no HIES. All eight TYK2-deficient patients displayed impaired but not abolished cellular responses to (a) IL-12 and IFN-α/ß, accounting for mycobacterial and viral infections, respectively; (b) IL-23, with normal proportions of circulating IL-17(+) T cells, accounting for their apparent lack of mucocutaneous candidiasis; and (c) IL-10, with no overt clinical consequences, including a lack of inflammatory bowel disease. Cellular responses to IL-21, IL-27, IFN-γ, IL-28/29 (IFN-λ), and leukemia inhibitory factor (LIF) were normal. The leukocytes and fibroblasts of all seven newly identified TYK2-deficient patients, unlike those of P1, responded normally to IL-6, possibly accounting for the lack of HIES in these patients. The expression of exogenous wild-type TYK2 or the silencing of endogenous TYK2 did not rescue IL-6 hyporesponsiveness, suggesting that this phenotype was not a consequence of the TYK2 genotype. The core clinical phenotype of TYK2 deficiency is mycobacterial and/or viral infections, caused by impaired responses to IL-12 and IFN-α/ß. Moreover, impaired IL-6 responses and HIES do not appear to be intrinsic features of TYK2 deficiency in humans.

The genetic control of immunity to Plasmodium infection.

BACKGROUND: Malaria remains a major worldwide public health problem with ~207 million cases and ~627,000 deaths per year, mainly affecting children under five years of age in Africa. Recent efforts at elaborating a genetic architecture of malaria have focused on severe malaria, leading to the identification of two new genes and confirmation of previously known variants in HBB, ABO and G6PD, by exploring the whole human genome in genome-wide association (GWA) studies. Molecular pathways controlling phenotypes representing effectiveness of host immunity, notably parasitemia and IgG levels, are of particular interest given the current lack of an efficacious vaccine and the need for new treatment options. RESULTS: We propose a global causal framework of malaria phenotypes implicating progression from the initial infection with Plasmodium spp. to the development of the infection through liver and blood-stage multiplication cycles (parasitemia as a quantitative trait), to clinical malaria attack, and finally to severe malaria. Genetic polymorphism may control any of these stages, such that preceding stages act as mediators of subsequent stages. A biomarker of humoral immunity, IgG levels, can also be integrated into the framework, potentially mediating the impact of polymorphism by limiting parasitemia levels. Current knowledge of the genetic basis of parasitemia levels and IgG levels is reviewed through key examples including the hemoglobinopathies, showing that the protective effect of HBB variants on malaria clinical phenotypes may partially be mediated through parasitemia and cytophilic IgG levels. Another example is the IgG receptor FcγRIIa, encoded by FCGR2A, such that H131 homozygotes displayed higher IgG2 levels and were protective against high parasitemia and onset of malaria symptoms as shown in a causal diagram. CONCLUSIONS: We thus underline the value of parasitemia and IgG levels as phenotypes in the understanding of the human genetic architecture of malaria, and the need for applying GWA approaches to these phenotypes.

Association study of genes controlling IL-12-dependent IFN-γ immunity: STAT4 alleles increase risk of pulmonary tuberculosis in Morocco.

BACKGROUND: Only a minority of individuals infected with Mycobacterium tuberculosis develop clinical tuberculosis. Genetic epidemiological evidence suggests that pulmonary tuberculosis has a strong human genetic component. Previous genetic findings in Mendelian predisposition to more severe mycobacterial infections, including by M. tuberculosis, underlined the importance of the interleukin 12 (IL-12)/interferon γ (IFN-γ) circuit in antimycobacterial immunity. METHODS: We conducted an association study in Morocco between pulmonary tuberculosis and a panel of single-nucleotide polymorphisms (SNPs) covering 14 core IL-12/IFN-γ circuit genes. The analyses were performed in a discovery family-based sample followed by replication in a case-control population. RESULTS: Out of 228 SNPs tested in the family-based sample, 6 STAT4 SNPs were associated with pulmonary tuberculosis (P = .0013-.01). We replicated the same direction of association for 1 cluster of 3 SNPs encompassing the promoter region of STAT4. In the combined sample, the association was stronger among younger subjects (pulmonary tuberculosis onset <25 years) with an odds ratio of developing pulmonary tuberculosis at rs897200 for GG vs AG/AA subjects of 1.47 (1.06-2.04). Previous functional experiments showed that the G allele of rs897200 was associated with lower STAT4 expression. CONCLUSIONS: Our present findings in a Moroccan population support an association of pulmonary tuberculosis with STAT4 promoter-region polymorphisms that may impact STAT4 expression.

CUBN and NEBL common variants in the chromosome 10p13 linkage region are associated with multibacillary leprosy in Vietnam.

Leprosy is caused by infection with Mycobacterium leprae and is classified clinically into paucibacillary (PB) or multibacillary (MB) subtypes based on the number of skin lesions and the bacillary index detected in skin smears. We previously identified a major PB susceptibility locus on chromosome region 10p13 in Vietnamese families by linkage analysis. In the current study, we conducted high-density association mapping of the 9.5 Mb linkage peak on chromosome region 10p13 covering 39 genes. Using leprosy per se and leprosy subtypes as phenotypes, we employed 294 nuclear families (303 leprosy cases, 63 % MB, 37 % PB) as a discovery sample and 192 nuclear families (192 cases, 55 % MB, 45 % PB) as a replication sample. Replicated significant association signals were revealed in the genes for cubilin (CUBN) and nebulette (NEBL). In the combined sample, the C allele (frequency 0.26) at CUBN SNP rs10904831 showed association [p = 1 × 10(-5); OR 0.52 (0.38-0.7)] with MB leprosy only. Likewise, allele T (frequency 0.42) at NEBL SNP rs11012461 showed association [p = 4.2 × 10(-5); OR 2.51 (1.6-4)] with MB leprosy only. These associations remained valid for the CUBN signal when taking into account the effective number of tests performed (type I error significance threshold = 2.4 × 10(-5)). We used the results of our analyses to propose a new model for the genetic control of polarization of clinical leprosy.

Age-dependent association between pulmonary tuberculosis and common TOX variants in the 8q12-13 linkage region.

Only a small fraction of individuals infected with Mycobacterium tuberculosis develop clinical tuberculosis (TB) in their lifetime. Genetic epidemiological evidence suggests a genetic determinism of pulmonary TB (PTB), but the molecular basis of genetic predisposition to PTB remains largely unknown. We used a positional-cloning approach to carry out ultrafine linkage-disequilibrium mapping of a previously identified susceptibility locus in chromosomal region 8q12-13 by genotyping 3,216 SNPs in a family-based Moroccan sample including 286 offspring with PTB. We observed 44 PTB-associated SNPs (p < 0.01), which were genotyped in an independent set of 317 cases and 650 controls from Morocco. A single signal, consisting of two correlated SNPs close to TOX, rs1568952 and rs2726600 (combined p = 1.1 × 10(-5) and 9.2 × 10(-5), respectively), was replicated. Stronger evidence of association was found in individuals who developed PTB before the age of 25 years (combined p for rs1568952 = 4.4 × 10(-8); odds ratio of PTB for AA versus AG/GG = 3.09 [1.99-4.78]). The association with rs2726600 (p = 0.04) was subsequently replicated in PTB-affected subjects under 25 years in a study of 243 nuclear families from Madagascar. Stronger evidence of replication in Madagascar was obtained for additional SNPs in strong linkage disequilibrium with the two initial SNPs (p = 0.003 for rs2726597), further confirming the signal. We thus identified around rs1568952 and rs2726600 a cluster of SNPs strongly associated with early-onset PTB in Morocco and Madagascar. SNP rs2726600 is located in a transcription-factor binding site in the 3' region of TOX, and further functional explorations will focus on CD4 T lymphocytes.

Crohn's disease susceptibility genes are associated with leprosy in the Vietnamese population.

A genomewide association study in Chinese patients with leprosy detected association signals in 16 single-nucleotide polymorphisms (SNPs) belonging to 6 loci, of which 4 are related to the NOD2 signaling pathway and are Crohn's disease susceptibility loci. Here, we studied these 16 SNPs as potential leprosy susceptibility factors in 474 Vietnamese leprosy simplex families. We replicated SNPs at HLA-DR-DQ, RIPK2, CCDC122-LACC1, and NOD2 as leprosy susceptibility factors in Vietnam. These results validated the striking overlap in the genetic control of Crohn's disease and leprosy.

Mycobacterial disease and impaired IFN-γ immunity in humans with inherited ISG15 deficiency.

ISG15 is an interferon (IFN)-α/ß-inducible, ubiquitin-like intracellular protein. Its conjugation to various proteins (ISGylation) contributes to antiviral immunity in mice. Here, we describe human patients with inherited ISG15 deficiency and mycobacterial, but not viral, diseases. The lack of intracellular ISG15 production and protein ISGylation was not associated with cellular susceptibility to any viruses that we tested, consistent with the lack of viral diseases in these patients. By contrast, the lack of mycobacterium-induced ISG15 secretion by leukocytes-granulocyte, in particular-reduced the production of IFN-γ by lymphocytes, including natural killer cells, probably accounting for the enhanced susceptibility to mycobacterial disease. This experiment of nature shows that human ISGylation is largely redundant for antiviral immunity, but that ISG15 plays an essential role as an IFN-γ-inducing secreted molecule for optimal antimycobacterial immunity.

Functional polymorphisms in IL13 are protective against high Schistosoma mansoni infection intensity in a Brazilian population.

BACKGROUND: IL-13 is a signature cytokine of the helper T cell type 2 (TH2) pathway which underlies host defense to helminthic infection and activates production of IgE in both parasitized populations and in urban settings after allergen exposure. METHODOLOGY/PRINCIPAL FINDINGS: Two functional polymorphisms in IL13, rs1800925 (or c.1-1111C>T) and rs20541 (or R130Q) were previously found to be associated with Schistosoma hematobium infection intensity. They have not been thoroughly explored in S. mansoni-endemic populations, however, and were selected along with 5 tagging SNPs for genotyping in 812 individuals in 318 nuclear families from a schistosomiasis-endemic area of Conde, Bahia, in Brazil. Regression models using GEE to account for family membership and family-based quantitative transmission disequilibrium tests (QTDT) were used to evaluate associations with total serum IgE (tIgE) levels and S. mansoni fecal egg counts adjusted for non-genetic covariates. We identified a protective effect for the T allele at rs20541 (P = 0.005) against high S. mansoni egg counts, corroborated by QTDT (P = 0.014). Our findings also suggested evidence for protective effects for the T allele at rs1800925 and A allele at rs2066960 after GEE analysis only (P = 0.050, 0.0002). CONCLUSIONS/SIGNIFICANCE: The two functional variants in IL13 are protective against high S. mansoni egg counts. These markers showed no evidence of association with tIgE levels, unlike tIgE levels previously studied in non-parasitized or atopic study populations.

A patient with tyrosine kinase 2 deficiency without hyper-IgE syndrome.

We describe a Turkish patient with tyrosine kinase 2 deficiency who suffered from disseminated Bacille Calmette-Guerin infection, neurobrucellosis, and cutaneous herpes zoster infection. Tyrosine kinase 2 deficiency should be considered in patients susceptible to herpes viruses and intramacrophage pathogens even in the absence of atopy, high serum IgE, and staphylococcal disease.

A mild form of SLC29A3 disorder: a frameshift deletion leads to the paradoxical translation of an otherwise noncoding mRNA splice variant.

We investigated two siblings with granulomatous histiocytosis prominent in the nasal area, mimicking rhinoscleroma and Rosai-Dorfman syndrome. Genome-wide linkage analysis and whole-exome sequencing identified a homozygous frameshift deletion in SLC29A3, which encodes human equilibrative nucleoside transporter-3 (hENT3). Germline mutations in SLC29A3 have been reported in rare patients with a wide range of overlapping clinical features and inherited disorders including H syndrome, pigmented hypertrichosis with insulin-dependent diabetes, and Faisalabad histiocytosis. With the exception of insulin-dependent diabetes and mild finger and toe contractures in one sibling, the two patients with nasal granulomatous histiocytosis studied here displayed none of the many SLC29A3-associated phenotypes. This mild clinical phenotype probably results from a remarkable genetic mechanism. The SLC29A3 frameshift deletion prevents the expression of the normally coding transcripts. It instead leads to the translation, expression, and function of an otherwise noncoding, out-of-frame mRNA splice variant lacking exon 3 that is eliminated by nonsense-mediated mRNA decay (NMD) in healthy individuals. The mutated isoform differs from the wild-type hENT3 by the modification of 20 residues in exon 2 and the removal of another 28 amino acids in exon 3, which include the second transmembrane domain. As a result, this new isoform displays some functional activity. This mechanism probably accounts for the narrow and mild clinical phenotype of the patients. This study highlights the 'rescue' role played by a normally noncoding mRNA splice variant of SLC29A3, uncovering a new mechanism by which frameshift mutations can be hypomorphic.

IL-12Rß1 deficiency in two of fifty children with severe tuberculosis from Iran, Morocco, and Turkey.

BACKGROUND AND OBJECTIVES: In the last decade, autosomal recessive IL-12Rß1 deficiency has been diagnosed in four children with severe tuberculosis from three unrelated families from Morocco, Spain, and Turkey, providing proof-of-principle that tuberculosis in otherwise healthy children may result from single-gene inborn errors of immunity. We aimed to estimate the fraction of children developing severe tuberculosis due to IL-12Rß1 deficiency in areas endemic for tuberculosis and where parental consanguinity is common. METHODS AND PRINCIPAL FINDINGS: We searched for IL12RB1 mutations in a series of 50 children from Iran, Morocco, and Turkey. All children had established severe pulmonary and/or disseminated tuberculosis requiring hospitalization and were otherwise normally resistant to weakly virulent BCG vaccines and environmental mycobacteria. In one child from Iran and another from Morocco, homozygosity for loss-of-function IL12RB1 alleles was documented, resulting in complete IL-12Rß1 deficiency. Despite the small sample studied, our findings suggest that IL-12Rß1 deficiency is not a very rare cause of pediatric tuberculosis in these countries, where it should be considered in selected children with severe disease. SIGNIFICANCE: This finding may have important medical implications, as recombinant IFN-γ is an effective treatment for mycobacterial infections in IL-12Rß1-deficient patients. It also provides additional support for the view that severe tuberculosis in childhood may result from a collection of single-gene inborn errors of immunity.

Accounting for genetic heterogeneity in homozygosity mapping: application to Mendelian susceptibility to mycobacterial disease.

INTRODUCTION: Genome-wide homozygosity mapping is a powerful method for locating rare recessive Mendelian mutations. However, statistical power decreases dramatically in the presence of genetic heterogeneity. METHODS: The authors applied an empirical approach to test for linkage accounting for genetic heterogeneity by calculating the sum of positive per-family multipoint LOD scores (S) across all positions, and obtaining corresponding empirical p values (EmpP) through permutations. RESULTS: The statistical power of the approach was found to be consistently higher than the classical heterogeneity LOD by simulations. Among 21 first-cousin matings with a single affected child, for five families linked to a locus of interest and 16 families to other loci, S/EmpP achieved a power of 40% versus 28% for heterogeneity LOD at an α level of 0.001. The mean size of peak linkage regions was markedly higher for true loci than false positive regions. The S/EmpP approach was applied to a sample of 17 consanguineous families with Mendelian susceptibility to mycobacterial disease, leading to the identification of two mutations in IL12RB1 and TYK2 from the largest of six linkage regions at p<10(-3). CONCLUSIONS: The S/EmpP approach is a flexible and powerful approach that can be applied to linkage analysis of families with suspected Mendelian disorders.

Germline CYBB mutations that selectively affect macrophages in kindreds with X-linked predisposition to tuberculous mycobacterial disease.

Germline mutations in CYBB, the human gene encoding the gp91(phox) subunit of the phagocyte NADPH oxidase, impair the respiratory burst of all types of phagocytes and result in X-linked chronic granulomatous disease (CGD). We report here two kindreds in which otherwise healthy male adults developed X-linked recessive Mendelian susceptibility to mycobacterial disease (MSMD) syndromes. These patients had previously unknown mutations in CYBB that resulted in an impaired respiratory burst in monocyte-derived macrophages but not in monocytes or granulocytes. The macrophage-specific functional consequences of the germline mutation resulted from cell-specific impairment in the assembly of the NADPH oxidase. This 'experiment of nature' indicates that CYBB is associated with MSMD and demonstrates that the respiratory burst in human macrophages is a crucial mechanism for protective immunity to tuberculous mycobacteria.